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Abstract 

It is shown that the 2x2 matrix Hamiltonian describing the dynamics of 
a charged spin 1/2 particle with g- factor 2 moving in an arbitrary, spatially 
dependent, magnetic field in two spatial dimensions can be written as the 
anticommuator of a nilpotent operator and its hermitian conjugate. Conse- 
quently, the Hamiltonians for the two different spin projections form partners 
of a supersymmetric quantum mechanical system. The resulting supersym- 
metry algebra can then be exploited to explicitly construct the exact zero 
energy ground state wavefunction for the system. Modulo this ground state, 
the remainder of the eigenstates and eigenvalues of the two partner Hamilto- 
nians form positive energy degenerate pairs. We also construct the spatially 
asymptotic form of the magnetic field which produces a finite magnetic flux 
and associated zero energy normalizable ground state wavefunction. 
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The calculation of the stationary states of an electron moving in an uniform magnetic 
field has been well studied for over 70 years and by now constitutes a classic textbook exam- 
ple The consequent Landau level energies and associated wavefunctions characterizing 
the dynamics in the two dimensional plane normal to the magnetic field have important 
applications in condensed matter systems ranging from the de Hass-van Alphen effect in 
metals || to the quantum Hall effect |§. Moreover, it has been observed {§- f| that when 
the particle has Lande factor g = 2, the system exhibits a supersymmetry with the Hamil- 
tonians corresponding to the spin up and spin down projections forming supersymmetric 
partners. The resulting supersymmetry algebra can be exploited to compute the exact zero 
energy ground state wavefunction as well as to show that the positive energy eigenvalues of 
the two Hamiltonians form degenerate pairs. It has also been previously noted || that the 
model describing the dynamics of a charged spin 1/2 particle with Lande factor g = 2 con- 
tinues to exhibit a supersymmetry even for arbitrary, spatially dependent, magnetic fields 
in two spatial dimensions. In the present note, we elaborate upon this observation and 
address some of its consequences. After explicitly constructing the supersymmetry charges 
and demonstrating the supersymmetry algebra, we use the construct to explicitly secure the 
exact form of the ground state zero energy wavefunction. The energy eigenvalues of the 
two supersymmetric partner Hamiltonians, which again correspond to the spin up and spin 
down projections, form degenerate pairs, save for the zero energy ground state of only one of 
the Hamiltonians. Finally, we discuss necessary conditions to produce normalizable ground 
states and finite magnetic field flux out of the plane. 

The Hamiltonian for a spin | particle of charge q and mass m moving in the x — y 
plane under the influence of an arbitrary, spatially dependent, static magnetic field in the 
^-direction, B = B(x, y)z, is given by 

(1 \ 

with the spin operator S = §03 = -| and the spin magnetic moment characterized 

Vo -1) 

by the Lande (/-factor. In two dimensions, the most general form for the magnetic vector 
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potential A is = 1,2) 



Ai(x,y) = diC(x,y) + e ij d j K{x,y) , 



(2) 



— * 

while the magnetic field is given by the curl of A as 



B(x,y) = etjdiAjfay) 



V 2 K(x,y) . 



(3) 



Adding a gradient to the magnetic vector potential leaves the magnetic field unaltered and 
corresponds to a gauge transformation. This implies that the scalar function C(x,y) is 
arbitrary while the prepotential K(x, y) is defined up to the addition of harmonic functions. 
Thus K(x, y) and K'(x, y) = K(x, y)+F(x+iy)+F*(x—iy) produce identical magnetic fields 
since the holomorphic functions F(x+iy) and F*(x—iy) automatically satisfy V 2 F(x+iy) = 
= V 2 F*(x — iy). As usual, all physical quantities such as the energies are necessarily gauge 
invariant. In the following, we choose to work in the Coulomb gauge defined by V • A = 
and with C = 0. The magnetic vector potential is then purely solenoidal and of the form 
Ai(x,y) = 6ijdjK(x,y), where K still retains the residual gauge freedom of the addition of 
F(x + iy) + F*(x - iy). 

It proves convenient to introduce complex coodinates x± = x ± iy, momenta p± = 
\ (Px T iPy) — f -£r^ = jd± and magnetic vector potential components A± = \ (A x ± iA y ) = 
^fidzpK, where the prepotential is now considered a function of x± so that K = K(x + ,x-). 
It follows that the magnetic field then can be written as B(x, y) = —Ad + d_K(x +) while 
the momenta conjugate to x± take the form 



Here = ^^-(H^ = are the respective Hamiltonians for the spin projections 




(4) 



and have the commutator [7r_,7r + ] = — 8-^d + d-K(x + , xJ) = 2-^B(x,y). 

When the Lande ^-factor takes the value g = 2, the Hamiltonian factorizes as 

/ iff \ l 1 7T + 7T_ \ 





+§( — §)■ This matrix Hamiltonian can further be written as the square of a hermitian 



operator Q as 
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1 „o „ 1 



— iir j 



H = -Q*;Q = — [ )=(Q)\ (6) 



2 xM 



Z7T_ 

An immediate consequence of this observation is that the energy spectrum is necessarily 
non-negative. 

/0 1\ /o 0' 

Using the Pauli matrices, cr + = |(<7i +icr 2 ) = ; a = \{(Ji — 1^2) 

\0 Oj \l 0, 

the operator Q can be further written as the sum Q = Q + Q* , where 

Q = —j=t^-0- ; Q f = —=-k + o+ (7) 



are two complex, nilpotent, Q 2 = = (<5^) 2 , supersymmetry charges. These charges, 
together with the Hamiltonian, obey the supersymmetry algebra |7|] [§] || 

{Q,Q}= ={Q\Q ] ) 
{Q,Q ] } = 2H 

[Q,H}= = [Q\H\. (8) 
The general normalizable eigenstate, ip^ n , has energy E^ n > and satisfies 

Hfifa n = ~^~^Tn = ^TnV'Tn ; > °- ( 9 ) 

Suppose if^ has the normalizable zero energy (E^q = 0) eigenstate ij)^ satisfying H^i/j^q = 
which implies that it-ip]® = 0. The vanishing commutator of Q and with H implies that 
7r_i7^ = H^tt_ and tt + H^ = H-^ti + . It follows that left multiplication of the eigenvalue 
equation by 7r_ then dictates that H^{ji_ip^ n ) = E^ n {ir _ip^ n ) , so that 7r_^ n , n > 0, is 
an ifj eigenstate with eigenvalue £"| n > 0. Note that the case n = does not give an 
eigenstate since it-ip]o = 0. Thus the normalizable ground state of is -?/>j.o = A ? |o 7r -V ; Ti> 
where iVj is a normalization constant, and has energy E^. Except for the zero energy 
eigenstate of iff, all the other eigenstates of iff and H± pair up with the same positive 
energy eigenvalues. This is a direct consequence of the supersymmetry. Thus if ip-\ n is an 
eigenstate of with eigenvalue E^ n , then %fj^ n = N± n 7r_^ n+ i , n > 0, is an eigenstate of 
Hi with eigenvalue = £q- n +i > 0. 



To explicitly construct the zero energy ground state of H^, we use that 

tt.^to = 2(-<9_ + i-d_K)^ = . (10) 

% c 

This differential equation is readily solved yielding 

^ (x + , x.) = N^U{x + )e^ K ^ x -) . (11) 

Here U (x + ) is an arbitrary complex function which reflects the degeneracy of the state. This 
zero energy eigenstate will be the system ground state provided it is normalizable. 
On the other hand, suppose the eigenvalue equation 

H^ ln = E in ^ in ; E ln >0 (12) 

admits the normalizable zero energy eigenstate, ^o, satisfying 7r + -?/^ = and E± = 0. 
Then an analogous argument gives that the normalizable ground state as ^-fo = ^TO 7r +'0J.i 
having positive energy E^ = E^. Once again, except for the zero energy eigenstate of 
all the other normalizable eigenstates of and pair up with degenerate positive energy 
eigenvalues. Thus ip^ n = N^ n n + i/ji n+ i is an eigenstate with eigenvalue E^ n = E± n+ i > 0. 
In this case, the zero energy ground state is gleaned from the condition 

tt+^o = 2(-«9 + - l ^d + K)iP l0 = (13) 
% c 

whose solution is 

lM* + , x.) = N i0 U(x-)e-& K <*+*-) (14) 

with Nio a normalization constant. This zero energy eigenstate will be the system ground 
state provided it is normalizable. If K is such that neither of the zero energy wavefunc- 
tions U(x + )e^ K ( x+ ' x ^ nor U(x-)e~^ K ^ x+ ' x ~^ are normalizable, then all states, including 
the ground states, have paired positive energies. This corresponds to a case of broken su- 
persymmetry. 
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These general results, obtained for an arbitrary x, y dependent magnetic field may also 
be applied to the example of a uniform magnetic field: B > 0. A corresponding prepotential 
is given by 

K(x + , x.) = ~\B( X -±^-r = -\ B y\ (15) 

Note that this choice of K corresponds to the asymmetric gauge choice (still in Coulonb 
gauge) A x = —By ; A y = 0. Recalling that K is only fixed up to additive functions 
F(x+) + F(x-), an alternate choice for K is K(x+.x-) = — \Bx + x^ = —\B(x 2 + y 2 ) 
which corresponds to the symmetric gauge A x = —\By ; A y = ^Bx. This is related to 
the asymmetric gauge choice by the holomorphic gauge transformation functions F(x+) = 
\Bx%. 

Working in the asymmetric (Landau) gauge, the state 

qB ( x -\-~ x — \ 2 qB 2 

iplo(x+,x_) = N^U(x + )e~^ [ ^~> = N^ U(x + iy)e~^ y (16) 

is the normalizable zero energy ground state. This is the familiar lowest Landau level 
eigenstate. Since any function can be expanded in terms of plane waves, we can choose 
U(x+) = e lkx+ and write 

ijMx,y)=N m e tkix+iy) e-^-y (17) 

where we have included a label k on the wavefunction which labels the degeneracy. Note 
that this wavefunction can be rewritten after completing the square as ipi k(x,y) = 
N-[oke lkx e~ %Kz( y ~ y °} 2 where y = —^0 and the overall normalization constant has been 
changed. For an unbounded system, the wavefunction, in this gauge, only exhibits a con- 
tinuum normalization in the x-direction, while the level is infinitely degenerate as there 
is no constraint on the allowed k and yo values. If the system is of finite extent, say a 
box of area L 2 , with L » ^|§, and periodic boundary conditions are imposed in the x 
direction so that ip]ok(x = —L/2,y) = ipiok(x = L/2,y), then the number of degenerate 
states is ^AA; = j^j^Ay — f-ip where we have used that Ay = L. Note that no par- 
ticular boundary condition has been imposed on the y direction. We are tacitly assuming 
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vanishing boundary conditions and although the wavefunction really does not vanish at 
y = ±|f, it is exponentially small there. On the other hand, the state e~ lkx ~ e^ y2 is clearly 
non-normalizable (in the continuum) so H± does not admit a zero energy normalizable eigen- 
state (all its eigenstates have positive energy). Note further that if the sign of B is reversed 
or if the sign of q is reversed (but not both) the role of and Hi are reversed and it is 
then which has the normalizable zero energy eigenvalue. 

In this example, one can, in fact, extract the entire spectrum. Noting that [^7^, ~^=\ — 
huj c , with u c = — , is precisely the commutation relation for the raising and lowering op- 
erators of a one dimensional simple harmonic oscillator, it follows that [ifj,7r + ] = huj c n + 
and consequently E^ n = nhu c ; n = 0,1,2,... (note absence of zero point energy), while 

n — iVj n (ir + ) n ^ o ; n — 0, 1, 2, ... Moreover, since and only differ by a constant, 
Hi — = —hu c , then Ei n = Ei n + Tiu c = n+i — (n + l)hu c ; n — 0, 1,2, ... and 

Note that in this uniform B field case, the spectrum can also be gleaned even in the case 
when g 7^ 2 and the supersymmetry is broken. For any value of g, the Hamiltonians Hi 
and Hi are simply modified by the additional constants —{g — 2)hu c and (g — 2)hu c . The 
resulting energy spectra for the spin up and spin down projections (assuming the spin up 
projection has the zero energy ground state when g — 2) are Ei n = hu c [n — (g — 2)] and 
Ei n = hu c [n+l+(g—2)]. Note that for g values other than 2, the ground state (n = 0) energy 
no longer vanishes relecting the broken supersymmetry. For the special value, g = 3/2, the 
two Hamiltonians are identical and their energy eigenvalues are ?Uv c (n+ |), while for g — 1, 
the roles of and Hi are reversed. Moreover, if g < 1 or g > 2, the spectrum contains a 
finite number of negative energy states. 

Let us return to the case of an arbitrary magnetic field and unbounded two dimensional 
space. In order for the zero energy eigenstate to be normalizable, it is necessary that the 
prepotential magnitude diverges at large distances. On the other hand, in order for the B 
field flux out of the x — y plane, 
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$= / dp p dipB= dtpipA^lp^ (18) 
Jo Jo Jo 

to be finite, the B field is required to fall off at large distances faster than \ = x 2+ y 2 ■ Here 
we have introduced the plane polar coordinates, p, (p, where x = pcos(p, y = psirup, so that 
x± = pe ±llfi , A v is the cp component of A and have employed Stokes' theorem to obtain the 
last equality. Taken together, in order to have both a normalizable zero energy ground state 
wavefunction and a finite B field flux, the prepotential must be proportional to powers of 
logarithms. A more detailed analysis reveals that the asymptotic behavior for K (p, (p) which 
meets both of the stated criteria can be written M as 



K( P , <p) ~ -^M^R = -^n(^) ; 1m > 0. (19) 



Here po is an arbitrary length scale in the logarithm. Since any value of po is as good as any 



other and thus no physical observable can depend on the po, we can simply take po = 



since it carries dimension length. The normalizable zero energy ground state wavefunctions 
are 

-^tn(Af) fl qq M . qq M 

^o(p,^)„~iV T onpV^e mc & = iV T0 „(— )^e m V-^, (20) 

1m 

where we have taken U(x + iy) — (x + iy) n = p n e m(p , with n a non-negative integer labeling 
the degeneracy. In fact, normalizability requires the parameters to satisfy n < — 1. Thus 
the ground state degeneracy, n + 1, is given by the largest integer less than 

Prepotentials depending on inp to a higher (lower) power give divergent (vanishing) B 
flux while producing normalizable (non-normalizable) zero energy wavefunctions. The above 
prepotential, in turn, corresponds to the asymptotic vector potential components: 

A p = A x cosip + AySintp ~ 

A<p = —A x sinip> + AyCosp ~ . (21) 

The resulting magnetic field strength then vanishes asymptotically 

5 (1) (p,V?) ~ (22) 
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while producing the magnetic flux 

$ = \q M - (23) 

A physical configuration exhibiting the above asymptotic vector potential is that of a 
vortex ||10|| . For a vortex characterized by the (quantized) mag neticflux , N = 1,2,..., 

it follows that = N and the ground state degeneracy is simply given by AT — 1. Note 
that this implies that the vortex solution carrying a single unit of magnetic flux, N — 1, 
does not produce a normalizable zero energy ground state in the unbounded plane. This is 
a straightforward consequence of the fact that the zero energy wavefunction in the presence 
of such a vortex falls asymptotically only as 1/p. An analogous result appears when the 
charged, spin 1/2 particle is described by a relativistic massless Dirac equation ||. In this 
case, the underlying connection between the integer characterizing the quantized magnetic 
flux of the vortex and the zero energy ground state degeneracy has a natural interpretation in 
terms of an index theorem denumerating the number of zero modes of the Dirac differential 
operator in a magnetic vortex background. Here we have explicitly seen that the zero energy 
ground state degeneracy is again related to the number of integer units of magnetic flux when 
the spin 1/2 particle is described by a non-relativistic Schrodinger equation where there is 
no index theorem connecting the two integer values. 

This work was supported in part by the U.S. Department of Energy under grant DE-FG02- 
91ER40681 (Task B). 
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